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Resumen 


En este trabajo se presentan los hallazgos obtenidos por medio de 
experimentos realizados en un canal de laboratorio, de sección 
transversal rectangular y paredes de vidrio, provisto de un fondo móvil 
constituido por un depósito de sedimentos uniformes no cohesivos, para 
el cual se simuló el inicio de movimiento. El objetivo perseguido 
consistió en estimar la rugosidad equivalente y las condiciones críticas 
para el inicio de movimiento. La rugosidad equivalente se obtuvo al 
contrastar los perfiles de velocidad obtenidos experimentalmente con el 
modelo teórico de Fuentes 8. Carrasquel (1978). La velocidad de corte 
se determinó a partir del esfuerzo turbulento de acuerdo con la 
distribución de los esfuerzos de Reynolds. La investigación incluyó el 
análisis de cinco diferentes tamaños de sedimento. Las velocidades 
instantáneas del flujo al centro del canal y sobre el fondo móvil fueron 
definidas por medio de la técnica óptica de velocimetría por rastreo de 
partículas (PTV). Lo anterior facilitó la estimación de la esperanza 
matemática puntual de la velocidad para intervalos del tirante de agua y 


con ello ajustar el perfil de velocidades de acuerdo con la ley logarítmica 
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de distribución de velocidad del flujo. Los resultados muestran que la 
rugosidad de fondo es mayor al diámetro medio de las partículas del 
lecho y que la razón de la rugosidad con respecto al tamaño medio del 
sedimento puede expresarse como una función del Reynolds de 
partículas. Se observó asimismo una disminución de la rugosidad 
equivalente justo después del inicio de movimiento de 26, 22 y 16% 
para lechos móviles con diámetro uniforme representativo (D) de 


partículas de 0.64, 0.93 y 1.68 mm, respectivamente. 


Palabras clave: sedimentos no cohesivos, rugosidad equivalente, inicio 


de movimiento, perfil de velocidades, PTV. 


Abstract 


The results obtained in experiments in a rectangular channel provided 
with glass walls and a movable bed in the center of the span are 
presented in this paper. Conditions for motion initiation of uniform non- 
cohesive sediment beds were obtained. The goal was to define the 
equivalent roughness and the critical conditions for initiation of motion. 
The bed roughness was obtained by comparing motion initiation velocity 
profiles with a theoretical velocity profile defined with the Fuentes 8 
Carrasquel (1978) model. The shear velocity was obtained with 
Reynolds stresses in a turbulent stress model. Five different size of 
sediments were used. Point velocities of flow at the channel and over 
the movable bed were obtained using the optical technique of Particle 
Tracking Velocimetry (PTV). From a cloud of measured velocities the 


average velocity profile was defined using the logarithmic law of flow 
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velocities. The results show that bed roughness is larger than the mean 
diameter of bed particles, and that the ratio of bed roughness over 
mean diameter can be expressed by means of the particle Reynolds 
number. It was observed a reduction of bed roughness just after motion 
initiation of 26, 22 and 16% for movable beds with uniform 


representative diameter (D) of 0.64, 0.93 and 1.68 mm, respectively. 


Keywords: Non-cohesive sediments, equivalent roughness, motion 


initiation, velocity profile, PTV. 


Recibido: 13/02/2019 
Aceptado: 08/07/2019 


Introducción 


La determinación de las condiciones de inicio de movimiento en lechos 
de sedimentos no cohesivos ha sido objeto de investigación desde Du 
Boys (1879). El trabajo más emblemático fue presentado por Shields 
(1936), quien recopiló datos experimentales y determinó la condición de 
inicio de movimiento con base en el esfuerzo cortante que actúa sobre el 
lecho. Posteriormente, Yalin 8 Karahan (1979) extienden la solución del 


inicio de movimiento para flujo laminar. Cabe destacar que la mayoría 
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de los resultados que se han obtenido para inicio de movimiento sugiere 
esfuerzos críticos que difieren en magnitud con los resultados de Shields 
(1936), o de Yalin € Karahan (1979), pero conservan la forma de las 
dos curvas propuestas respectivamente (Buffinton, 1999; Cao et al., 
2006; Chien 8 Wan, 1998; García-Flores 8 Maza-Álvarez 1997; Lick et 
al., 2004; Ling, 1995; Smith, 2004). Investigaciones complementarias 
han sido documentadas sobre la influencia de las condiciones de 
empaque de las partículas del lecho y de las diferencias de tamaño de 
éstas con respecto al inicio de movimiento, apoyándose en conceptos 
estadísticos para explicar la variabilidad del fenómeno flujo-fondo móvil 
(Dancey et al. 2002; Garde et al., 2006; Papanicolaou, 2002). Por otro 
lado, diversos autores (Hunziker 8 Jaeggi, 2002; Kuhnle, 1993; 
Sarmiento 8 Falcon, 2006) subrayan que los lechos de ríos están 
constituidos por lo general por mezclas de arena o granos de diferentes 
tamaños, de tal forma que de modo prioritario debe construirse un 
modelo matemático que explique el problema para fondo con partículas 
de tamaño uniforme y relacionar, o extrapolar, dicha solución para 


lechos formados por partículas de granulometría extendida. 


La condición para inicio de movimiento en un fondo móvil se 
presenta por lo regular en términos del esfuerzo cortante, o de la 
velocidad media, teniendo ventajas y desventajas cada uno de estos 
enfoques. El enfoque de velocidad media del flujo para explicar el inicio 
de movimiento ha sido adoptado en la práctica de la ingeniería por su 
sencillez. Así, modelos como los propuestos por Levy (1956) o 
Goncharov (1962) difieren sobre todo en el tamaño de partícula 


representativa de la rugosidad del lecho en su respectivo modelo. Sin 
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embargo, autores como Chien 8 Wan (1998), y García-Flores 8 Maza- 
Álvarez (1997) señalan, con base en el análisis del efecto de la 
turbulencia, que la velocidad media del flujo es un criterio menos 
eficiente que el del esfuerzo cortante para explicar de manera correcta 
el fenómeno del inicio de movimiento. Además, Marsh et al.(2004) 
concluyen que el empleo de la variable velocidad media de flujo como la 
principal variable explicativa del fenómeno de inicio de movimiento del 
fondo debe ser analizada o validada con respecto a los esfuerzos de 
corte generados por ésta. En consecuencia, estos últimos autores 
proponen el uso de la variable potencia del flujo como una forma 
relativamente nueva de modelación del inicio de movimiento del lecho 
(Marsh et al., 2004). 


Cabe recordar que la rugosidad equivalente es el parámetro que 
representa la altura de las asperezas del lecho móvil. Nikuradse (1933) 
propone que la rugosidad equivalente del fondo es igual al tamaño 
promedio de las partículas para lechos de arenas uniformes. Sin 
embargo, resultados posteriores han demostrado que la rugosidad 
equivalente es mayor al tamaño medio de las partículas que lo 
componen (Van-Rijn, 1982). Además, para fondo no uniforme 
(granulometría extendida), la rugosidad equivalente del fondo debe ser 
representada por alguno de los diámetros más grandes de la mezcla de 
partículas. En la literatura especializada, la rugosidad equivalente es 
sugerida como un diámetro representativo de las partículas del lecho 
ponderado por un factor mayor que la unidad. El problema por resolver 
consiste, entonces, en determinar la magnitud de la rugosidad 


equivalente del fondo móvil y definir el diámetro de partícula 
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representativo del lecho, justificando el criterio utilizado (Chien € Wan, 
1998; Maza-Álvarez 8 García-Flores, 1996). Por otro lado, en la práctica 
actual de la ingeniería, continúan utilizándose las fórmulas de Manning y 
Chézy para la solución de problemas de resistencia al flujo. En efecto, 
los coeficientes de rugosidad utilizados en tales modelos se obtienen 
como una función del tamaño representativo de las asperezas de la 
pared sobre la cual ocurre el flujo. En consecuencia, las correctas 
magnitudes de la rugosidad equivalente del fondo y de la resistencia al 
flujo para explicar el inicio de movimiento del lecho son de suma 
importancia (Ramesh et a/., 2000; Wu € Wang-Sam, 1999). 


La interacción entre la pared y fondo del canal, o sección mojada 
de río, con el flujo se refleja en la distribución de velocidades. Se ha 
demostrado que la mejor aproximación corresponde con un perfil 
logarítmico de velocidades en la vertical, que es representado por un 
modelo matemático que, a su vez, es una función de la rugosidad y del 
flujo presentes, donde el esfuerzo cortante está implícito en la 
distribución logarítmica (Maza-Álvarez 8 García-Flores, 1984). De esta 
forma, el perfil de velocidades logarítmico permite definir una velocidad 
para el inicio de movimiento del fondo relacionada con el esfuerzo 
cortante y la rugosidad equivalente. Por ello, el objetivo perseguido en 
este trabajo consistió en medir perfiles de velocidad a partir de la 
técnica óptica de velocimetría por rastreo de partículas (PTV) (Salinas- 
Tapia, 2007) y obtener combinaciones de esfuerzo cortante-velocidad 
media generadoras del inicio de movimiento del fondo móvil, para luego 
contrastarlas con perfiles teóricos de velocidad de acuerdo con el 


modelo de Fuentes 8 Carrasquel (1978). De esta forma ha sido posible 
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deducir valores de la rugosidad equivalente de fondo para cinco 
diferentes tamaños de arena y relacionarlos con la condición de inicio de 


movimiento del lecho móvil. 


En la actualidad, gracias a las técnicas ópticas (no intrusivas), 
como PTV, es posible la medición de las fluctuaciones de la velocidad o 
bien de la distribución del esfuerzo de Reynolds. Lo anterior representa 
una alternativa para estimar el esfuerzo cortante en situaciones donde 
la medición de la pendiente del fondo es incierta, como el caso de los 
ríos. La velocidad de corte y, por ende, el esfuerzo cortante total, se 
obtiene con las fluctuaciones de la velocidad promedio en la altura de 
flujo (Yang 8: McCorquodale, 2004; Lemmin € Rolland, 1997). 


La distribución, o perfil, de velocidades promedio es resultado de 
las interacciones del flujo con las fronteras que lo contienen. Un modelo 
obtenido con base en el análisis teórico, ampliamente sustentado, que 
relaciona parámetros con significado físico es la ley logarítmica de 
distribución de velocidades de Prantl-von Karman (Anthony et al., 


1998), el cual está representado por la Ecuación (1): 
¿InqZ+ (1) 
Ux k Ks Pp 


Donde u es la velocidad promedio; ux, la velocidad de corte (m/s); 
k, la constante de von Karman (k = 0.41 para flujo con baja carga 
suspendida); y, la altura de flujo (m); v, la viscosidad cinemática 


(m/s), y p es una función de la naturaleza hidráulica de la pared o 
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función de rugosidad. La constante de von Karman (k) representa la 
relación proporcional entre la magnitud de la turbulencia del flujo y la 
variación de la velocidad promedio. La velocidad de corte (ux) es un 
parámetro que representa al esfuerzo cortante del fondo o promedio To 
en unidades de velocidad (u* = (ro/p)*?, p = densidad del fluido). La 
función de rugosidad O es igual con la diferencia entre una constante de 
integración KB y el logaritmo natural del parámetro conocido como 
Reynolds de rugosidad (Rex); esto es: 0 = KB -— In(Rex), donde Rex = 
Uuxk./v (para ks = altura de las asperezas de la pared del fondo). Nótese 
que al sustituir Hd = KB - In(uxks/v) en la Ecuación (1), la ley logarítmica 


de distribución de velocidades toma la forma de la Ecuación (2): 
“12 
E In ra +6 (2) 


La constante KB es función de la naturaleza hidráulica del fondo, 
expresada en términos de Rex. Fuentes 8. Carrasquel (1978) proponen la 


Ecuación (3) para obtener la magnitud de B. 


M0 (3) 


La altura que representa la rugosidad equivalente de las asperezas 
del fondo (ks) continúa siendo tema de interés en investigaciones 


recientes. En efecto, en canales abiertos con fondo plano constituido por 
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partículas móviles y no cohesivas (arenas y/o gravas) se ha 
documentado ampliamente que ks es igual al tamaño representativo de 
los granos del fondo (D;¡) multiplicado por una constante de 
proporcionalidad a (donde a > 1). Ejemplo de lo anterior son los 
trabajos de Ackers € White (1973), quienes proponen un valor de ks = 
1.25 D35; Engelund 8 Hansen (1967) sugieren un ks = 2 Des; O bien 
Einstein, en 1950 (Einstein, 1950), quien sugirió un ks = Des. Van-Rijn 
(1982) hace una recopilación de datos experimentales de distintos 
autores y concluye que la rugosidad de un fondo plano y móvil varía de 
la 10 veces el Doo, proponiendo como rugosidad promedio a ser usada 
Ks = 3 Doy. Señala también que la rugosidad no depende de la condición 
de transporte en el fondo cuando éste permanece plano. Así, de acuerdo 
con los trabajos precedentes, una solución al problema planteado en 
este trabajo compara los resultados con los métodos que usan la 
constante de proporcionalidad (a) y el tamaño representativo de los 
granos del lecho móvil (D;). La condición de inicio de movimiento ha sido 
adoptada como criterio para definir el tamaño representativo de un 
fondo compuesto por una mezcla de tamaños como normalmente se 
presenta en condiciones naturales (ríos y canales). Para el caso de fondo 
uniforme, el tamaño promedio de partículas será considerado como 
adecuado para representar el tamaño de las partículas del lecho. Aquí 
cabe resaltar que en la modelación hidráulica, por lo regular se utilizan 
coeficientes de rugosidad como n (Manning) y C (Chézy), los cuales 


pueden estimarse como una función de ks (Chien 8 Wan, 1998). 


Con respecto al inicio de movimiento de un fondo móvil, la 


condición crítica para la cual los granos empiezan a ser arrastrados por 
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la corriente ha sido establecida como una relación entre dos parámetros 
adimensionales: el parámetro de Shields (0) y el Reynolds de partículas 
(Rep) (Re, = uxD/v). El parámetro de Shields (0) representa la relación 
entre la acción de la corriente (esfuerzo cortante en el fondo) y el peso 
sumergido de los granos. Se tiene como valor de 8 = To/D(ys-y), donde 
Ys es el peso específico del grano; y es el peso específico del fluido 
(agua). Cabe diferenciar entre los parámetros Rex« y Rep, dado que el 
primero se refiere a la naturaleza hidráulica de la pared del fondo, 
mientras el segundo considera la acción de la corriente sobre los granos 
del fondo. Considerando que ks > D, o bien que a > 1, se tiene entonces 


que normalmente Rex> Re),. 


El conocido diagrama de Shields (1936) es la representación 
gráfica de la condición crítica para inicio de movimiento del lecho. Así, 
en una extensión del trabajo de Yalin 8 Karahan (1979) presentan un 
diagrama análogo que incluye la solución del inicio de movimiento del 
lecho para condiciones de flujo laminar. Por otro lado, García-Flores 8 
Maza-Álvarez (1997) proponen ecuaciones en función de R€p, que 
ajustan a los diagramas propuestos por Shields, y Yalin y Karahan, 
respectivamente. Es de resaltar que en la región correspondiente a la 
transición de pared lisa a pared hidráulicamente rugosa, los diagramas 
tanto de Shields como de Yalin y Karahan presentan la forma 
coloquialmente conocida como silla de montar. Según Shields, la 
relación entre la dimensión de la subcapa viscosa y el tamaño de los 
granos es tal que para granos contenidos dentro del espesor de la 
subcapa viscosa se requiere una fuerza de arrastre mayor para romper 


el equilibrio que guardan éstos. La condición crítica, o de inicio de 
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movimiento del lecho es un parámetro considerado en varios métodos 
para estimar el transporte de fondo. Por ejemplo, en el modelo 
propuesto por Meyer-Peter 8 Muller (1948), la magnitud del parámetro 
9 = 0.045 es adoptada cuando Re, > 60. 


Si bien el concepto de esfuerzo cortante crítico es más adecuado 
bajo el análisis del equilibrio individual de los granos del fondo, en la 
práctica la velocidad media del flujo resulta ser un criterio de mayor 
simplicidad para estimar el inicio del movimiento o arrastre del fondo. 
Modelos clásicos para calcular la velocidad media del flujo, que ocasiona 
el inicio de movimiento del lecho, consideran tanto la distribución de 
velocidades en la vertical como la condición crítica de estabilidad de las 
partículas del fondo en términos del parámetro de Shields (8). Uno de 
estos modelos clásicos corresponde al propuesto por Chien 8 Wan 


(1998) y está representado por la Ecuación (4). 
UL xR 
——= =5B.75 [f(Rep) log(12.27)7- (4) 
Ys Yop Ñ 
J y 


Donde U. es la velocidad media crítica para inicio de movimiento 
del lecho (m/s); R, el radio hidráulico de la sección transversal (m); Kks, 
el valor de rugosidad equivalente de las partículas (m); g, el valor de la 
aceleración de la gravedad (m“/s); D, el diámetro representativo de las 
partículas del lecho móvil (m); x, la constante en función de la relación 


entre en espesor de la subcapa viscosa Ó (m) y la rugosidad k;, y f(Rep) 
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(N/m?) es el parámetro de Shields (8) en la condición de inicio de 
movimiento del lecho. Con base en los resultados de Shields (1936) 
para Re, mayores a 60, el valor de 68 es considerado como constante y 
cercano a 0.045. Así, la Ecuación (4) se simplifica, y el modelo de Chien 


8, Wan (1998) puede ser representado por la Ecuación (5): 


Uc 


=1.219 [log(12.27)£] (5) 


Ys Y 
y 92 


Por su parte, Einstein (1950) propone la determinación de 
existencia de una pared rugosa cuando k./S > 10 y x = 1. Igualmente 
señala que si ks/ 0 < 0.25, entonces se está en presencia de una pared 
lisa y x toma el valor de x = 0.3Rex. El valor máximo posible de x (x = 
1.6) ocurre cuando se está en zona de transición y psa para cuando k; / 
9 es igual a 1.0. Es importante señalar que varios modelos encontrados 
en la literatura especializada son de la forma de la Ecuación (5). Entre 
éstos se encuentran los modelos de Goncharov (1962) y de Levy 
(1956), representados por la Ecuación (6) y la Ecuación (7), 


respectivamente: 


8.8h 


“e 1.06 log) (6) 


Ys Y 
y 22 
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Donde h es el tirante del flujo (m), Doo y Dos (m) son, 
respectivamente, el diámetro que representa que 90 o 95% de las 


partículas del fondo móvil es menor a estos valores. 


==1,4 log (7) 


La principal diferencia entre estos modelos es la definición de la 
rugosidad de los granos (k;). Los citados autores suponen que esa 
rugosidad está representada por el tamaño mayor de partículas en el 
lecho, es decir, Doo O Dos, respectivamente (Chien 8 Wan, 1998; García- 
Flores 8 Maza-Álvarez, 1997). Las tres expresiones anteriores incluyen 
siempre la relación entre la profundidad (h), y/o el radio hidráulico (R), 
así como la rugosidad de los granos (ks). Recordando que por definición 
el esfuerzo cortante promedio en el lecho está representado por Tp = 
yRs, donde s es la pendiente de fricción considerada. La relación R/ks 
tiene en cuenta la acción de la corriente, es decir, la potencia del 
escurrimiento (R = h en canales donde B >> h). Para estimar el 
esfuerzo cortante sobre el fondo, en condiciones donde la pendiente de 
fricción no puede conocerse directamente, se utilizan otros métodos. En 
ríos, por ejemplo, la velocidad de corte es deducida al validar la 
aplicación de la ley logarítmica de distribución de velocidades. Otra 
forma de deducir la velocidad de corte, e indirectamente el esfuerzo 
cortante promedio, es validar la distribución de los esfuerzos de 


Reynolds 7, con el empleo del modelo representado por la Ecuación (8): 
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7,=pu? (1- 2) (8) 


Por definición, T, = pu'v”, donde u” es la variabilidad de la 
velocidad con respecto al promedio en dirección del flujo; v”, la 
variabilidad de la velocidad promedio en la dirección vertical; el valor de 
y se encuentra entre 0 y h. La distribución anterior ha sido validada para 
flujo turbulento en canales abiertos. Al respecto, Lemmin € Rolland 
(1997) compararon el esfuerzo cortante total calculado con base en la 
pendiente de fricción (s) y el obtenido directamente al extrapolar las 
distribuciones del esfuerzo de Reynolds hasta el fondo, y ampliaron la 
aplicación para flujo en cauces naturales (Best et al., 1997; Song € 
Graf, 1996). Por último, es importante señalar que la técnica PTV es no 
intrusiva y permite medir las fluctuaciones de la velocidad en diferentes 
puntos del tirante del flujo para establecer, con base en dichas 
mediciones, la velocidad de corte ux (Adrian, 1991; Salinas-Tapia et al., 
2006). 


Metodología 
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Diagrama procedimental 


El diagrama de la Figura 1 describe el proceso metodológico que se llevó 


a Cabo en este trabajo. 
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Modelación del inicio 
de movimiento en 
canales con fondo 

móvil 


Método propuesto por 


Modelos clásicos - 
el presente trabajo 


Medición del perfil de 
velocidad para inicio 
de movimiento por 
PTV 


Esfuerzo cortante Velocidad crítica 
diagrama de Shields media, Levy (1956), 
(1936) Goncharev (1962) 


Definición empírica de Medición de velocidad 
rugosidad k. = f(a)D; de corte por PTV (u”v') 


Formulación del perfil 
teórico de velocidad 
(Fuentes €: Carrasquel 
(1978) 


Deducción de k; por 
contraste de perfiles 
de velocidad teóricos y 
observados 


Figura 1. Diagrama del proceso metodológico empleado en la 


investigación. 
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Montaje experimental 


Las dimensiones del canal de laboratorio donde se realizaron las 
simulaciones del inicio de movimiento del lecho de fondo móvil son 1.20 
m de longitud, 0.10 m de ancho y altura de las paredes laterales de 
0.15 m. A 0.60 m de la entrada del canal se ubica un depósito de fondo 
móvil de 0.025 m de profundidad y 0.25 m de longitud. Las partículas 
de arena, o fondo móvil, se colocaron en dicho depósito verificando que 
el lecho presentara siempre una configuración plana al inicio de cada 
experimento. En la salida del canal, una compuerta permitió controlar la 
velocidad del flujo, manteniendo un caudal constante durante cada 
experimento. En la Figura 2 se muestra el esquema del dispositivo 
canal-depósito de fondo móvil y la configuración del sistema de 


velocimetría por rastreo de partículas empleado (PTV). 
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Láser 


pulsado Lente 


D 
Región de 
Análisis Perfil de 
> “welocidad 


Hoja de luz 


1.20 m. su Fondo 
2 2. móvil 
Canal y Región de análisis 


Láser _ Espejo 
pulsado Lente 
Hoja de luz 


Cámara CCD 


' h Fondo 
—— móvil 


Procesador 0.10 m. 
Sección Transversal 


Figura 2. Configuración del montaje experimental de iluminación, 
captura y procesamiento de imágenes en la región de flujo sobre fondo 


móvil. 
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Como fondo móvil se utilizó arena de mina, cuya densidad y peso 


específico son ps = 2 543 kg/m* y ys = 24 950 N/m?, La arena fue 


tamizada para obtener cinco diferentes tamaños uniformes (lechos 


móviles) D. La Tabla 1 contiene el resumen de las aberturas mínima y 


máxima (Dmn Y Dmx) del tamiz utilizado para obtener los cinco tamaños 


D. El promedio de Dimn Y Dmx se considera como el diámetro uniforme en 


cada caso. Se muestra, igualmente, la clasificación correspondiente del 


material de cada lecho con respecto al sistema de clasificación del 


American Geophysical Union (AGU). 


Tabla 1. Fondo móvil y clasificación AGU (American Geophysical Union) 


del sedimento que lo conforma. 


Tamaño Tamaño 
Tamaño o 
a mínimo máximo Clasificación AGU del 
diámetro 
Fondo móvil ] según según sedimento del fondo 
uniforme: D 
tamiz Dm. | tamiz Dmx móvil 
(m) 
(m) (m) 
s1 0.00038 0.00033 0.00043 Arena media 
S2 0.00064 0.00043 0.00085 Arena gruesa 
S3 0.00093 0.00085 0.00100 Arena gruesa 
S4 0.00168 0.00100 0.00236 Arena muy gruesa 
S5 0.00318 0.00236 0.00400 Grava fina 
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En cada simulación efectuada del inicio de movimiento del lecho 
de fondo, la velocidad del flujo se incrementó de modo gradual con 
(del 


aproximadamente en cada paso). Después de calibrar las condiciones 


aberturas mínimas de la compuerta orden de 0.5 mm 
que generan el inicio de movimiento del lecho móvil se efectuaron los 
experimentos a partir de velocidades donde si bien son próximas a la 
condición crítica, se mantiene en reposo de la totalidad de los granos del 
lecho. Cada experimento fue extendido hasta condiciones evidentes de 
transporte de fondo. El inicio de movimiento se consideró por criterio 
visual cuando para una condición de flujo dada es incuestionable que 
Cada 


experimento incluyó dos tipos de resultados correspondientes a las 


hay un movimiento generalizado y constante de granos. 


condiciones “antes de movimiento” y “después de movimiento”. La Tabla 
2 muestra un resumen de los parámetros hidráulicos medidos en cada 
experimento. Dicha tabla contiene las alturas de flujo h, las velocidades 
medidas al principio y al final de cada etapa durante la simulación, así 
como sus correspondientes Res (número de Reynolds del fluido) y Fr 


(número de Foude). 


Tabla 2. Parámetros hidráulicos al principio y final de las simulaciones 


de inicio de movimiento del lecho de fondo. 


Fondo móvil: ] Altura | Número de | Número de 
Velocidad Caudal 
tamaño, Etapa del . de Reynolds | Froude Fr . 
, . media U, . . promedio 
clasificación | experimento flujo h del flujo = de 
(m/s) 12 | 9 (m”/s) 
AGU (m) |Res¿ = Uh/v| U/(gh) 
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s1: D= 
0.00038 m, 


arena media 


Inicio, antes 
de 


movimiento 


0.199 


0.049 


8089 


0.29 


Final, 
después de 


movimiento 


0.298 


0.035 


8739 


0.51 


0.00103 


Inicio, antes 
de 


movimiento 


0.247 


0.075 


15300 


0.29 


Final, 
después de 


movimiento 


0.276 


0.064 


14599 


0.35 


0.00180 


S2:D= 
0.00064 m, 


arena gruesa 


Inicio, antes 
de 


movimiento 


0.256 


0.042 


8852 


0.40 


Final, 
después de 


movimiento 


0.320 


0.033 


8862 


0.56 


0.00105 


Inicio, antes 
de 


movimiento 


0.287 


0.061 


14431 


0.37 


Final, 
después de 


movimiento 


0.354 


0.051 


15027 


0.50 


0.00178 


S3: D= 
0.00093 m, 


arena gruesa 


Inicio, antes 
de 


movimiento 


0.234 


0.042 


8216 


0.36 


0.00100 
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Final, 
después de 


movimiento 


0.302 


0.034 


8432 


0.53 


Inicio, antes 
de 


movimiento 


0.287 


0.063 


14876 


0.37 


Final, 
después de 


movimiento 


0.379 


0.045 


14178 


0.57 


0.00175 


S4: D = 
0.00168 m, 


Inicio, antes 
de 


movimiento 


0.304 


0.055 


13996 


0.41 


arena muy 


gruesa 


Final, 
después de 


movimiento 


0.441 


0.041 


14828 


0.70 


0.00174 


Inicio, antes 
de 


movimiento 


0.298 


0.036 


8998 


0.50 


S5: D = 3.18 


mm, grava 


Final, antes 
de 


movimiento 


0.350 


0.031 


8899 


0.64 


0.00108 


fina 


Inicio, antes 
de 


movimiento 


0.491 


0.035 


14377 


0.84 


Final, antes 
de 


movimiento 


0.693 


0.027 


15485 


1.35 


0.00180 
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Como se observa en la Tabla 2, los experimentos se hicieron para 
un rango de velocidades medias de 0.199 a 0.693 m/s; y relaciones 
ancho (B)/tirante (Ah) de 1.3 a 3.7. Se destaca que la relación B/h 
durante los experimentos es baja con respecto a otros estudios 
documentados en la literatura; este hecho se discute en el análisis de 
los resultados obtenidos. Con base en la Ecuación (9), propuesta por 
Kirkgóz 8 Ardichoglu (1997), para definir la longitud L necesaria de un 
canal y garantizar que la región turbulenta incluyera toda la profundidad 


del canal, se verificó en todos los experimentos. 


=76-0.0001 (9) 


El perfil de velocidades se dedujo al ajustar un modelo de 
regresión logarítmico con la información de las velocidades promedio u y 
sus correspondientes profundidades y. Las esperanzas matemáticas de 
cada coordenada (u, y) fueron obtenidas para intervalos de 1 mm de 
altura de flujo. Por lo tanto, u, y son los promedios de las velocidades 
instantáneas u; y las correspondientes alturas de flujo y; contenidas en el 
intervalo vertical de amplitud 0.001 m. El plano, o lámina, de 
iluminación del sistema óptico se ubica en cada medición sobre el centro 
del canal y del fondo móvil, formando una región de análisis de 0.08 m 
de longitud en la dirección del flujo y altura igual con el tirante h. Un 
total de 100 imágenes doble pulsadas fueron capturadas a razón de 60 
imágenes por segundo. El procesamiento de las imágenes facilitó, por 
medio de la técnica PTV, la medición de un promedio de 7 000 vectores 
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de velocidad (u;, vi) durante cada repetición del experimento para cada 
par de coordenadas (x;, y;) de la ventana del canal analizada y para cada 
paso de abertura de la compuerta en cada una de las distintas 


simulaciones. 


Resultados 


Perfiles de velocidad 


En la Figura 3 se tienen ejemplos del perfil de velocidades obtenido para 
un caso o medición por cada fondo móvil utilizado. El gráfico mostrado 
en la esquina inferior derecha de la Figura 3 incluye el comparatrivo de 
los perfiles de velocidad correspondientes a los cinco ejemplos 


presentados en la misma figura. 
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Figura 3. Velocidades instantáneas obtenidas por la técnica PTV y 
ajustes de sus respectivos promedios con la ley logarítmica de 
distribución de velocidades; a, b, c, d, e) un ejemplo por cada fondo 


móvil utilizado; f) comparativo de los perfiles de velocidad mostrados. 
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Como se observa en los ejemplos de la Figura 3, el perfil de 
velocidades se obtuvo al aplicar a los promedios u, y el modelo de 
regresión logarítmico para el cual se obtiene una curva de la forma 


expresada por la Ecuación (10): 


y=Aebou (10) 


Donde y es la altura de flujo (m); u, la velocidad promedio para la 
altura y (m/s); A y bo son las constantes del modelo de regresión. En 
cada caso, la región tomada en cuenta para aplicar el modelo de 
regresión fue acotada con base en la ventana donde se obtiene el mejor 
coeficiente de determinación R* posible, definiéndose así la región de 
flujo en la vertical donde se garantiza que la distribución de velocidades 
sigue la ley logarítmica (región de flujo turbulento). La amplitud de 
dicha región depende sobre todo de la profundidad h en cada medición, 
profundidad que, como se mencionó, fue controlada por pequeños 
incrementos de apertura de la compuerta a la salida del canal. Una serie 
experimental o simulación de inicio de movimiento implicó un conjunto 
de mediciones desde una condición de flujo próxima al inicio de 
movimiento del lecho hasta una condición donde el arrastre de los 
granos del fondo fue indiscutible. Así, en cada simulación se tienen 
mediciones “antes de movimiento” y “después de movimiento”. Los 
gráficos insertos en la Figura 4 muestran ejemplos de una de las 
simulaciones por cada fondo móvil, excepto para S5, para el cual se 


presentan las tres simulaciones realizadas. 
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Figura 4. Perfiles de velocidad durante la simulación de inicio de 
movimiento. Para S1, S2, S3, S4 se muestra una simulación por fondo 


móvil; para S5 se observan las tres simulaciones realizadas. 
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Esfuerzo cortante y rugosidad del fondo 


El parámetro ux puede obtenerse a partir de la constante bo de la curva 
de regresión obtenida (Ecuación (10)) para las velocidades medidas 
experimentalmente, es decir ux*x = k/bo. Lo anterior se demuestra al 
comparar la Ecuación (2) y Ecuación (10). La constante A corresponde a 
la altura de flujo para la cual teóricamente, de acuerdo con el modelo 
logarítmico, se tiene u = O. Este parámetro es conocido como rugosidad 
hidráulica y comúnmente designado con el símbolo yo. La velocidad de 
corte ux considerada para estimar el esfuerzo cortante To se tomó con 
base en la distribución de los esfuerzos de Reynolds (Ecuación (8)). Por 
lo tanto, Uux es estimada con base en la Ecuación (11) (Chien 8 Wan, 
1998): 


a 
ux=(u'v”) (2= A (11) 


Donde u% es el promedio de la variabilidad de la velocidad en la 
dirección del flujo; v”, el promedio de la variabilidad de la velocidad en 
la dirección vertical; ¡, el tamaño del intervalo considerado, y M es el 


número total de intervalos o dovelas en los que se divide la profundidad 
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h. Las respectivas variabilidades promedio se obtienen a partir de la 


Ecuación (12) y la Ecuación (13). 


u=- de (u,-u)” (12) 


i 2 
v="Y4 (v;-v) (13) 


Donde u, y v; son las velocidades puntuales (de los intervalos 
considerados) en la dirección del flujo y la vertical, respectivamente; z 
es el número de datos o vectores de velocidad contenidos en el intervalo 
analizado. La velocidad de corte (ux) obtenida con respecto a los 
esfuerzos de Reynolds representa esfuerzos de corte que tienen mayor 
congruencia con las simulaciones efectuadas. Para las simulaciones con 
el fondo móvil S5 (D = 0.0032 m), el esfuerzo cortante logrado al 
considerar ux con respecto a los esfuerzos de Reynolds es correcto, dado 
el hecho de no haber observado inicio de movimiento del lecho. Por el 
contrario, si el esfuerzo cortante se tomara con respecto a ux del perfil 
logarítmico, se obtienen magnitudes mucho mayores que el valor 
correspondiente al esfuerzo crítico según las curvas de Shields, y Yalin y 
Karahan. La Figura 5 muestra el contraste entre la velocidad de corte 
obtenida por la ley logarítmica de velocidades ux(log) y la velocidad de 


corte por medio de los esfuerzos de Reynolds ux(rs). 
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a) antes de movimiento o 


0.04 qe > 0.04 


b) después de movimiento 


O|." 


= E > 
un un $ 
=> e 

£ 0.03 £, 0.03 o ES 
a ne 24 MÚÓ xS2 
— AO 
£ £ AA A S3 
3 3 a e 

0.02 0.02 e 054 

0.01 a + A + z + S + A + E 1 0.01 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.01 0.02 0.03 0.04 


ux(log): [m/s] ux(log): [m/s] 


Figura 5. Comparación: velocidad de corte obtenida por la ley 
logarítmica ux(log) versus velocidad de corte obtenida por los esfuerzos 
de Reynolds ux(rs); a) antes de inicio de movimiento del lecho; b) 


después de movimiento del lecho. 


Con base en la transformación logarítmica de la Ecuación (2) y la 
Ecuación (10), y  despejando el término /n(y) se obtienen, 
respectivamente, la Ecuación (14) y la Ecuación (15). Las relaciones bo 
= k/ux* e In(A) = In(ks) - Bk son evidentes al contrastar la Ecuación (14) 


y la Ecuación (15): 


In(y)=b,u+In(A4) (14) 


In(y)= Lu+In(ks)-Bk (15) 


Tecnología y ciencias del agua, ISSN 2007-2422, 11(2), 105-157. DOI: 10.24850/j-tyca-2020-02-03 
135 


2020, Instituto Mexicano de Tecnología del Agua 
Tecnología y “ Open Access bajo la licencia CC BY-NC-SA 4.0 
Ciencias¿Agua 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 
Así, con base en el modelo de Fuentes y Carrasquel (Ecuación 
(3)), la constante B es obtenida como una función de Rex (Rex = uxky/v). 
Por lo tanto, la magnitud de ks para la cual se cumple la identidad /n(A) 


= In(ks) - Bk se obtiene al resolver la Ecuación (16): 


In)=In +0 La (16) 


La Ecuación (16) resulta de sustituir el modelo de Fuentes y 
Carrasquel en /n(A) = In(ks) - Bk. La constante A es obtenida con base 
en el ajuste de las velocidades promedio medidas en el modelo 
logarítmico (Ecuación (10)); y ux se consideró conforme con los 
esfuerzos de Reynolds al calcular el promedio del producto de los 
cuadrados medios de las fluctuaciones de la velocidad. Así, al resolver la 
Ecuación (16) se obtuvo la magnitud de la rugosidad equivalente del 


fondo móvil (k;). 


Rugosidad equivalente del fondo móvil antes y después 


del movimiento 
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Tres condiciones del experimento afectan la amplitud de la región de 
flujo turbulento: la distancia desde la entrada del canal hasta la lámina 
de iluminación (0.60 m); la relación B/h, y el remanso producido por la 
compuerta. Por tales razones, la región (o ventana) donde la ley 
logarítmica describe de forma adecuada la distribución de velocidades 
promedio fue definida en cada caso con respecto al mejor coeficiente de 
determinación encontrado (R?). La región de aplicación de la ley 
logarítmica se incrementó conforme mayor fue la apertura de 
compuerta. En otras palabras, la amplitud de la altura de flujo en la que 
la ley logarítmica es significativa, pues la ventana de análisis aumentó 
con el tirante durante las simulaciones. En la Tabla 3 se muestra el 
resultado para cada fondo móvil simulado, donde se incluye el promedio 
de los coeficientes de determinación (R?) de los experimentos 
correspondientes. 


Tabla 3. Región de flujo donde la ley logarítmica describe de manera 


adecuada la distribución de velocidades en la vertical del flujo. 


Región que 
corresponde al R? promedio en los 
Fondo móvil - tamaño mejor coeficiente | perfiles de velocidad 
de determinación observados 
R? 
S1 - D = 0.00038 m [0.05 - 0.65]h 0.93 
S2 - D = 0.00064 m [0.04 - 0.83]h 0.94 
S3 - D = 0.00093 m [0.04 - 0.76]h 0.95 
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S4 - D = 0.00168 m [0.06 - 0.88]h 0.97 


S5 - D = 0.00318 m [0.08 - 1.00]h 0.98 


En los promedios presentados en la Tabla 3 y en la determinación 
de la rugosidad de fondo se consideraron sólo perfiles de velocidad para 
R? > 0.90. Para las simulaciones con fondo S1, se obtuvo la menor 
región de correspondencia con la ley logarítmica y los menores 
coeficientes de determinación; mientras que para el fondo S5 se tienen 
la mejor correlación y la mayor amplitud de la región de flujo que 
corresponde a ley logarítmica. La amplitud y correlación con la ley 
logarítmica dependieron de la profundidad del tirante durante el 
experimento. De acuerdo con la Ecuación (9), la mínima longitud 
necesaria para el desarrollo completo del flujo turbulento es L = 2 m, 
longitud mayor que la distancia física experimental entre la entrada del 
canal e inicio del fondo móvil. Sin embargo, para las simulaciones con 
fondo S5, la amplitud de la región donde la ley logarítmica es adecuada, 
incluye toda la profundidad, lo que podría considerarse como un flujo 


turbulento completamente desarrollado. 


Magnitudes de rugosidad efectiva (k./D) de 0.91 a 5.47 y de 1.14 
a 5.46 se obtuvieron para las condiciones “antes de movimiento” y 
“después de movimiento”, respectivamente. En la Figura 6 se presentan 
los gráficos de (k-/D) con respecto al Reynolds de partícula (Rep) y al 
parámetro de Shields (0). De acuerdo con los resultados obtenidos, la 
rugosidad equivalente de fondo (k«/D) disminuye conforme aumentan 


ambos parámetros. 
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Figura 6. Rugosidad equivalente (ks/D) vs. Rep: a) antes de 
movimiento; b) después de movimiento, rugosidad equivalente (k./D) 


vs. 8; c) antes de movimiento; d) después de movimiento. 


Los gráficos de la Figura 6 incluyen todos los resultados k./D. Los 
promedios por cada serie experimental o simulación se muestran en la 
Figura 7. Esta figura incluye también la curva del modelo de regresión 
del tipo ks/D = K*Rep” (donde K y n son los coeficientes del modelo de 


regresión). En la figura 7 se observan los promedios k«/D con respecto a 
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8 y dos diferentes modelos de regresión: recta, ky/D = MO + b (donde M 
y b son los coeficientes del modelo de regresión), y curva, ky/D = K208” 
(para cuando K2 y m son los coeficientes del modelo de regresión). Las 
curvas de regresión observadas en la Figura 6 omiten los resultados 
obtenidos para el fondo Si. Las rugosidades equivalentes de fondo 
(k¿/D) obtenidas para Si no corresponden al comportamiento oO 
tendencia mostrada por los resultados obtenidos para los otros cuatro 
fondos (S2-S5), lo cual es fácilmente observado en los gráficos 
mostrados de la Figura 5 y Figura 6. Los menores coeficientes de 
determinación (R?*) y las menores amplitudes de ventana de análisis, 
donde la ley logarítmica es adecuada, corresponden a las simulaciones 
S5, lo cual es congruente con la discrepancia de k./D obtenida para el 


mismo lecho con respecto al resto de éstos. 
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Figura 7. Modelos de regresión obtenidos para el promedio de la 
rugosidad equivalente (k.“/D) antes y después del inicio de movimiento: 
a) ks/D vs. Rep; b) ks/D vs. 6. 


El promedio ponderado de la rugosidad equivalente de fondo 
(ks/D), y sus respectivos máximos y mínimos se muestran en la Tabla 4. 
Con respecto a los promedios antes y después de movimiento, se 
observa que hay una reducción de ks/D para los fondos S2, S3 y S4; en 


tanto que para S5, el resultado ks/D es prácticamente el mismo antes y 
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después de movimiento. Sin embargo, como ya se discutió, para el 
fondo S5, las consideraciones del análisis realizado no son estrictamente 
válidas (menores correlaciones y regiones de aplicación de la ley 
logarítmica). Para S2, S3, S4, la reducción de la rugosidad equivalente 
“después de movimiento” fue de 26, 22 y 16%, respectivamente, 
comparada con la rugosidad equivalente en condición “antes de 
movimiento”. Es importante señalar que en investigaciones que tratan el 
tema de la rugosidad en condiciones de transporte de fondo se concluye 
que el arrastre de los granos aumenta la rugosidad equivalente del lecho 
(Van-Rijn, 1982; Camenen et al., 2008). Sin embargo, se subraya que 
ese aumento de rugosidad equivalente del lecho corresponde a una 
situación sostenida de arrastre de fondo y que en el caso de las 
simulaciones realizadas en el presente trabajo k“y/D “después de 
movimiento” se refiere a un evento inmediato al inicio de movimiento 


del lecho. 


Tabla 4. Rugosidad equivalente de lecho (k«/D) antes y después de 


movimiento (valores promedio y rango de los resultados). 


Kks/D: antes de k./D: después de 


movimiento movimiento 


Iced haced ica mázca 


Fondo móvil - 


tamaño 


S1 - D= 0.00038 m 4.12-5.47 4.25-5.56 
S2 - D= 0.00064 m 3.16 2.42-3.96 2.34 1.60-3.23 
S3 - D= 0.00093 m 2.65 1.73-3.39 2.08 1.65-2.49 
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S4 - D = 0.00168 m 2.02 1.60-2.34 1.70 1.14 - 2.16 


S5 - D= 0.00318 m 1.33 0.90-1.93 MA 


Discusión 


En la Figura 8 (gráfico superior) se muestran los resultados obtenidos en 
términos de Re, vs. 8, donde se incluyen ambas condiciones de la 
simulación (antes y después de movimiento). En la misma figura 
(gráfico inferior) se muestra el resultado correspondiente a la condición 
para la cual se observó el inicio de movimiento en cada simulación. Para 
S5 no se produjo inicio de movimiento; los resultados mostrados en la 
gráfica inferior de la Figura 8 corresponden al máximo esfuerzo cortante 
o parámetro de Shields (8) producido en las tres simulaciones 


realizadas. 
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Ta) simulaciones: antes y después de movimeinto A Sniclds 
= = ="Yalin 8 Karahan" 

D Si (antes de movimiento) 
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x  S2 (antes de movimiento) 
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movimiento) 

1 10 100 a  S3 (después de 
Rep: uxD/v movimiento) 


0.01 


b) condición crítica ——- "Shields" 

= += "Yalin 8 Karahan" 
B” si 
x  S2 


a  Ss3 
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Figura 8. Comparación de los resultados obtenidos con respecto a las 
curvas de Shields, y Yalin y Karahan: a) resultados antes y después de 


movimiento; b) inicio de movimiento observado en los experimentos. 


Esfuerzos cortantes críticos y por tanto parámetros de Shields 
mayores a los propuestos por las curvas de Shields, y de Yalin y 
Karahan fueron obtenidos para S1, S2 y S3; mientras que para S4, los 
resultados corresponden a esfuerzos de corte menores al teórico 
calculado según ambas curvas. De acuerdo con el gráfico inferior de la 
Figura 7, de haberse producido experimentalmente el inicio de 
movimiento para el fondo S5, la condición crítica se ubicaría por encima 
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de las curvas referidas, generándose al unir las condiciones críticas una 


tendencia similar a la de dichas curvas. 


Como velocidad media crítica se consideró la correspondiente al 
último perfil de velocidades, determinado justo antes del inicio de 
movimiento en cada simulación. En la Tabla 5 se muestran los 
resultados experimentales (PTV) con respecto a la condición de inicio de 
movimiento, y las magnitudes de los parámetros críticos teóricos 
calculados por el criterio de Shields (esfuerzo cortante) y los modelos de 
Chien y Wan, Goncharov, Levy, así como el de Fuentes y Carrasquel, 
respectivamente. En los tres modelos de velocidad media crítica se 
utilizó ks = D para los cinco lechos móviles analizados. Mientras que 
para el modelo de Fuentes y Carrasquel, se consideró ks igual con el 
resultado observado en la Tabla 4 (antes de movimiento) para S2, S3, 
S4 y S5, pero para Si se tomó ks = 8.2Rep-0.38. El parámetro de 
Shields (0) se consideró con respecto a la magnitud teórica según las 
expresiones de Maza -Alvarez 8 Garcia-Flores (1984). En el caso del 
fondo S5 (D = 3.18 mm,), se incluyen la velocidad media y el esfuerzo 
cortante máximos medidos experimentalmente, los cuales, como se 
indicó de manera oportuna, no corresponden al inicio de movimiento 


sino a condiciones sin movimiento. 


Tabla 5. Resultados de PTV y comparación con otros autores. 


Fondo móvil Parámetros r. Velocidades medias - inicio de 


- tamaño críticos criterio movimiento: U, (m/s) 


(m) experimentales de k. =D k. (PTV) 
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(PTV) Shields 
(N/m?) | Chien 
U. Te Fuentes 8 
8 Goncharov | Levy 
m/s |(N/m?) Carrasquel 
Wan 

s1:D= 

0.234 0.33 0.21 0.235 0.241 0.245 0.242 
0.00038 
S2:D= 

0.285 0.44 0.31 0.263 0.282 0.288 0.280 
0.00064 
S3:-D = 

0.309 0.60 0.47 0.313 0.338 0.336 0.348 
0.00093 
S4:D= 

0.362 0.69 1.02 0.404 0.400 0.400 0.460 
0.00168 
S5:D= 

0.702 1.94 2.36 0.521 0.470 0.472 0.632 
0.00318 


*** PTY = máximos medidos experimentalmente sin movimiento observado 


Las velocidades medias críticas experimentales (obtenidas por 
PTV) para S1 y S2 son cercanas a las calculadas por los cuatro modelos 
considerados. No obstante, para estos mismos fondos, los esfuerzos 
cortantes críticos experimentales son mayores a los obtenidos por el 
modelo de Shields. En los casos S3 y S4, las velocidades críticas 
calculadas son mayores que las observadas de modo experimental, 
mientras los correspondientes esfuerzos cortantes experimentales son 
mayores en el caso S3 y menores para S4. En el lecho S5, aún cuando 
experimentalmente la máxima velocidad medida es mayor que la 


calculada por cualquiera de los cuatro modelos considerados, el esfuerzo 
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cortante que corresponde a esta velocidad máxima (0.703 m/s) es 
menor que el propuesto por el modelo de Shields. También para esta 
condición se tiene Fr > 1 (ver Tabla 1), en tal combinación velocidad- 
profundidad, la energía específica del flujo durante el experimento 
disminuye conforme aumenta la velocidad y decrece la profundidad. El 
gráfico mostrado en la Figura 9 contiene el comparativo entre los 
parámetros críticos experimentales y los criterios teóricos considerados 


contenidos en la Tabla 5. 


Inicio de movimiento: experimental vs teóricos 


2.5 
mel = e PTV 
2.0 , 
oChiens8Wan 
Ti.5 x Goncharov 
E 
LA 
g A Levy 
e 1.0 
DFuentes8gCarrasquel 
0.5 mS5 (sin movimiento) 
0.0 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
U. [m/s] 


Figura 9. Velocidad media y esfuerzo de corte críticos vs. criterio de 
Shields (modelos de velocidad crítica: Chien y Wan, Goncharov, Levy, y 


Fuentes y Carrasquel). 
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Conclusiones 


De acuerdo con los hallazgos encontrados en este trabajo para 
simulaciones de inicio de movimiento de un fondo móvil generado por 
un flujo en canal, se tiene que la rugosidad equivalente para un fondo 
plano, constituido por partículas uniformes no cohesivas, debe 
considerarse un factor de ponderación mayor a uno, que afecte al 
tamaño medio de los granos. Además, dicho valor de ponderación es 
variable con respecto a la clasificación del tamaño de las partículas del 
fondo móvil, obteniéndose un mayor valor para arenas medias que para 
aquellas clasificadas como arenas gruesas y gravas finas. La rugosidad 
equivalente de fondo móvil se analizó con relación al Reynolds de 
partículas (Rep) y al parámetro de Shields (8), obteniéndose, para 
ambos casos, altos coeficientes de determinación entre los modelos 
empíricos y los resultados experimentales. Al considerar la diferencia 
entre las condiciones “antes de movimiento” y “después de 
movimiento”, se concluye que la respuesta del lecho es, en términos 
hidráulicos, correspondiente con una pared de menor rugosidad justo 
después de ser rebasado el estado límite de inicio de movimiento. Se 
subraya que tal hallazgo responde exclusivamente a situaciones 
inmediatas después del inicio de movimiento donde el fondo permanece 


plano. 
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La magnitud de la rugosidad equivalente, considerada al aplicar 
criterios de velocidad media crítica, influye en el cálculo de ésta. En el 
caso de arenas gruesas y gravas finas, se observaron las mayores 
diferencias entre las velocidades medias críticas calculadas y las 
medidas experimentalmente con ayuda de la técnica PTV. Para los 
fondos utilizados en las simulaciones de tales clasificaciones 
granulométricas, se obtuvo una sobreestimación de la velocidad media 
crítica cuando se consideró como rugosidad equivalente de fondo al 


tamaño medio de los granos. 


La técnica PTV permitió obtener cantidades considerables de 
mediciones de velocidad instantánea, a partir de las cuales se pudieron 
calcular los esfuerzos de Reynolds y perfiles de velocidad promedio, 
estimándose de manera indirecta el esfuerzo cortante en el lecho móvil 
del canal. El esfuerzo de corte obtenido, al acotar la región de flujo, con 
base en la ley logarítmica, es comparable con el logrado a partir de los 
esfuerzos de Reynolds en condiciones de flujo turbulento y subcrítico. 
Para régimen supercrítico, el esfuerzo cortante obtenido por los 
esfuerzos de Reynolds resultó ser un parámetro más congruente con los 
resultados arrojados con respecto a la condición de no movimiento 
observado para la grava fina. Cabe señalar que tanto la distancia desde 
la entrada del canal al depósito de fondo móvil y la longitud física del 
canal afectaron de forma directa la amplitud de la región turbulenta en 


la zona de interés. 


Al comparar los resultados de esta investigación con los de otros 
investigadores, se puede recomendar para calcular la rugosidad efectiva 
el uso de una constante de proporcionalidad (a) entre 1.33 y 4.92, 
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afectando el diámetro medio del sedimento (Dso) de acuerdo con la 
Tabla 4. 


Con base en la discusión antes expuesta, simulaciones similares 
para un canal de mayor longitud son recomendables en futuros 
experimentos. Sin embargo, los esfuerzos críticos obtenidos 
corresponden con el comportamiento típico de las curvas de inicio de 
movimiento clásicas. Simulaciones como las aquí presentadas pueden 
realizarse bajo las mismas consideraciones aplicadas, con objeto de 
inferir la magnitud de la constante de proporcionalidad (a) y el tamaño 
representativo de los granos del fondo (D), a partir del criterio de inicio 
de movimiento. Por último, la realización de una serie de experimentos 
con mezclas de tamaño de partículas no cohesivas (granulometría 
extendida) serían ¡igualmente recomendables para extrapolar la 


aplicabilidad de los hallazgos encontrados. 


Nomenclatura 


f: Constante o función de rugosidad en el perfil logarítmico de 


velocidades. 


bo, A: Constantes del modelo de regresión del perfil de velocidades 


experimental. 
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D: Tamaño de partícula representativo (mm). 
D35: Tamaño de partícula mayor al 35% (mm). 
Ds5o: Tamaño de partícula mayor al 50% (mm). 
Des: Tamaño de partícula mayor al 65% (mm). 
Doo: Tamaño de partícula mayor al 90% (mm). 
Do5: Tamaño de partícula mayor al 95% (mm). 
Dmn: Tamaño mínimo de partículas (mm) (máxima abertura de malla). 
Dmx: Tamaño máximo de partícula (mm) (mínima abertura de malla). 
g: Aceleración de la gravedad (9.81 m/s?). 
h: Tirante del escurrimiento (cm). 
k: Constante de von Karman (0.41). 
ks: Rugosidad de fondo (mm). 
R: Radio hidráulico (mm). 
Rx: Reynolds de rugosidad (adimensional). 
Rep: Reynolds de partículas (adimensional). 
s: Pendiente de fricción (m/m). 
S1, S2, S3, S4, S5: Fondos móviles de diámetro uniforme. 
u: Velocidad media puntual (cm/s). 
Uux. Velocidad de corte (cm/s). 
u”“: Variabilidad de la velocidad u (cm/s). 


U¿: Velocidad media crítica (cm/s). 
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U,: Velocidad media flujo (cm/s). 

v”: Fluctuación de la velocidad v (cm/s). 

V.: Velocidad media en la dirección vertical (cm/s). 
y: Altura de flujo (m). 

X: Parámetro o función de rugosidad de Einstein (adimensional). 
vs: Peso específico de las partículas (N/m?). 

v: Viscosidad cinemática del agua (m/s). 

9: Parámetro de Shields (adimensional). 

ps: Densidad de las partículas (kg/m?). 

To: Esfuerzo cortante en el fondo (N/m?). 

T.: Esfuerzo cortante crítico (N/m?). 


T: Esfuerzo de Reynolds o esfuerzo cortante (N/m?). 
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